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CALCULATION  OF  THE  DOSE  RATE  EMANATING 
FROM  GAMMA  RADIATION  OF  A  GAS  JET 
PROPAGATING  THROUGH  THE  SURFACE  LAYER  OF 
THE  ATMOSPHERE 

0.  V.  Zakharov,  V.  V.  Kovalenko,  and 
V.  M.  Kolobashkin. 

Using  the  principle  of  the  semiempirical  theory  of  turbulent  dif¬ 
fusion,  work  [1]  has  solved  the  problem  of  a  radioactive  mixture  prop¬ 
agating  in  the  surface  layer  of  the  atmosphere. 

Let  us  examine  the  problem  of  calculating  the  dose  rate  on  the 
earth's  surface  which  emanates  from  a  gas  jet,  whose  concentration  of 
activity  is  distributed  in  accordance  with  the  solution  derived  in 
work  [1],  for  the  various  energies  of  the  gamma  quanta  and  varied 
meteorological  conditions. 


Fig.  1.  Geometry  for  calculating  the  dose  rate. 


Dose  rate  P(x0)  at  point  xQ  (Fig.  1)  from  a  volume  source  can  be 
calculated  by  the  formula 

where  q(x,  y,  z)  is  the  activity  concentration  obtained  in  work  [1], 
eurie/mJ;  -  differential  gamma  constant  of  an  isotope,  p-m  /(h* 
curie);  r  -  distance  from  the  volume  element  dV  (m^)  to  the  point  xQ, 
m;  -  linear  coefficient  of  attenuation  of  gamma  radiation  with  the 
energy  Eit  m“*;  Bp(p^r)  -  dose  accumulation  factor  proposed  in  work 
[2]  for  E^iO.5  MeV  and  in  work  [3]  for  E^O.S  MeV;  the  sum  is  taken 
over  all  the  gamma  lines  of  a  given  isotope. 

Integration  is  performed  with  respect  to  the  entire  jet  volume  V, 
where  the  activity  concentration  is  different  from  zero.  Fig.  1  shows 
the  relationship  between  the  rectangular  and  spherical  coordinates. 

If  points  xQ,  in  the  dose  rates  are  being  calculated,  lie  at  the 
projection  of  the  jet's  axis  on  the  earth's  surface,  then,  due  to  the 
symmetry  of  the  problem  relative  to  the  XOZ  plane,  the  $  integral  is 
selected  in  the  range  from  0  to  it  and  it  is  doubled. 


The  function  P^(x0,  S^),  which  numerically  equal  to  the  dose  rate 
in  p/h  created  by  the  gamma  radiation  of  the  jet  flowing  from  a  source 
with  power  Q0=l  curie/s  with  the  wing  velocity  v=l  m/s,  was  calculated 
for  a  long-lived  isotope  with  curie  •  Thus,  below  we  will  refer 

to  the  function  P^CXq,  E^)  as  the  reduced  dose  rate. 

The  upper  limit  of  integration  with  respect  to  r  -  R=800  m  (using 
the  computer,  we  have  determined  that  the  contribution  when  integrating 
for  R28OO  m  yields  an  addition  of  less  than  IK  of  P^  when  R=800  m). 

The  absence  of  the  jet  at  points  x£0  was  taken  into  account  by  the  fact 
that  the  integrand  was  assumed  to  be  zero  for  points  xQ<800  m,  for 
x=xQ+r  sin  0  cos  $s0. 

Integral  (1)  was  calculated  by  a  numerical  method  using  the 
Simpson's  formula.  In  our  case,  the  initial  subdivision  steps  are 

and  =20  m. 

After  calculating  with  the  initial  step,  this  step  is  reduced  one-half 

and  the  integral  is  calculated  again,  etc.  The  calculations  using 

each  variable  end  when  the  next  decrease  of  the  subdivision  step  by 

one-half  yields  the  relative  correction  of  the  integral  value  of  less 
_ 2 

than  10  .  The  reduced  dose  rate  was  calculated  on  a  digital  computer 

for  the  energies  of  the  gamma  quanta  E^  (MeV}=0. 10-2.5  and  six  values 
of  the  atmosphere-stability  parameter  e=-0.3;  -0.2;  -0.1;  0.0;  0.1;  and 
0.2. 

The  reduced  dose  rates  calculated  thus  for  the  most  frequently  en¬ 
countered  indifferent  stratification  of  the  atmosphere  (e=0)  are  rep¬ 
resented  in  Fig.  2  by  dashed  curves  as  functions  of  distance  to  the 
source  for  gamma  radiation  with  E^sl  MeV.  In  the  same  figure,  the 
solid  lines  depict  the  reduced  surface  concentration  q(xQ)  as  a  func¬ 
tion  of  distance  to  the  source,  which  is  numerically  equal  to  the  con¬ 
centration  of  the  long-lived  gaseous  isotope  in  curie/m^  from  a  source 
with  power  Q0=l  curie/s  at  the  wind  velocity  v=l  m/s  at  the  points 
lying  at  the  projection  of  the  jet's  axis  on  the  earth's  surface.  The 
reduced  surface  concentration  has  the  dimensionality  [m~2]. 

The  values  of  the  reduced  dose  rate  are  shown  on  the  right  axis  of 
the  ordinates,  the  values  of  the  reduced  surface  concentration  -  on 
the  left  axis  of  the  ordinates. 

The  reduced  dose  rate  P^  has  a  slight  dependence  on  the  energy  of 


gamma  radiation:  in  the  energy  range  from  0.1  to  2.5  MeV  it  changes 
by  ±30%  from  the  values  with  E  i  =  1  MeV.  The  dependence  of  on  the 
distance  to  the  source  is  the  same  for  all  energies  of  the  range  ex¬ 
amined. 

The  following  formula  can  be  used  to  calculate  the  dose  rate  P  from 
a  source  with  power  Q  (curie/s)  for  a  specific  isotope  with  the  wind 
velocity  v  (m/s): 

f  •  (2} 

Taking  into  account  the  slight  dependence  of  P^Xq,  Bj)  on  energy  Bj 
with  an  error  j£30J,  the  expression  for  the  dose  rate  from  the  gamma 
radiation  of  a  specific  isotope  can  be  written  as 

'w  — f  ww.f.  <3; 

o  m2 

where  is  the  complete  gamma  constant  of  the  isotope,  ^  curie* 


The  effect  of  the  radioactive  decay  of  the  isotope  with  A/0  can  be 
taken  into  account  by  multiplying  the  value  of  P(xQ),  calculated  by 
formula  (2)  or  (3),  by  •  .  Calculation  on  the  computer  shows  that 
in  taking  into  account  the  decay  of  the  isotope 

(ji*. 

*<*.» 1  ' 

where  P(xQ,  A)  is  the  dose  rate  at  distance  xQ  from. the  source  for  an 
isotope  with  the  decay  constant  A;  P(xQ)  -  analogous  dose  rate  from  a 
long-lived  isotope  -  yields  diffference  from  the  values  of  P(xQ,  A) 
calculated  by  formula  (1)  with  the  addition  of  the  term,  which  takes 
into  account  a  decay  of  not  over  8%  for  the  isotopes  with  a  half-decay 
period  T ^22:20  min  and  xQ>l-10^  m  with  the  wind  velocity  vsl  m/s.  This 
difference  is  less  for  high  wind  velocities. 

It  can  be  seen  from  Pig.  2  that  the  position  of  the  maximum  of  sur¬ 
face  concentration  of  the  isotope  does  not  coincide  with  the  position 
of  the  maximum  of  the  surface  dose  rate  from  the  gamma  radiation  of  the 
same  isotope.  This  noncoincidence  of  the  maxima  is  observed  at  all  the 
energies  of  the  gamma  quanta  and  meteorological  parameters  observed. 

The  maximum  of  the  surface  dose  rate,  in  the  case  of  a  raised  source 
(hrfO),  is  not  clearly  expressed  and  is  located  at  short  distances  from 


the  source.  We  have  the  same  thing  for  the  stable  (e>0)  and  unstable 
(e<0)  temperature  stratifications  of  the  atmosphere. 

In  all  these  cases ,  the  maxima  of  surface  concentration  of  gases 
are  located  considerably  further  from  the  source. 

These  relationships  between  the  maxima  of  surface  dose  rate  and 
surface  concentration  of  gases  must  be  taken  into  account  when  deter¬ 
mining  the  dimensions  of  the  safety  zones  for  the  concerns  in  atomic 
industry. 
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